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ABSTRACT 



Rece nt models o f varia tions of the Sun's structure with the 11-year activity 
cycle by ISofia et al.l (120051 ) predict strong non-homologous changes of the radius 
of subsurface layers, due to subsurface magnetic fields and field-modulated tur- 
bulence. According to their best model the changes of the surface radius may be 
1000 times larger than those at the depth of 5 Mm. We use f-mode oscillation 
frequency data from the MDI instrument of Solar and Heliospheric Observatory 
(SOHO) and measurements of the solar surface radius variations from SOHO and 
ground-based observatories during solar cycle 23 (1996-2005) to put constraints 
on the radius changes. The results show that the above model overestimates the 
change of the radius at the surface relative to the change at 5Mm. 



Subject headings: Sun: helioseismology — Sun: oscillations — Sun: activity - 
Sun: interior 
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Introduction 



For the last 30 years, the question of whether the solar radius changes with time, and 
with the 11-year activity cycle, in particular, has been the subject of many controversies. 
Using different techniques, several series of measurements have been made leading so far to 
inconsistent results. Even if a strolabe measurements can be con s idered as extracted from 
the same statisti c al pop ulation ( Badache-Damiani fc Rozelotll2006 ). Laclare et al. ( 1996 ) and 
Reis Neto et al.l (120031) reported a variation of the solar radius in antiphase with the solar 



cycle while iNoell (12004 ) showed the opposite. In addition, drift-time measurements made 
photoelectrically from CCD trans its at Izana and Locarno d id not show cycle- dependent 
variations in exces s of ±50 mas3 ( Wittmann fc Bianda 2000 ). Using SOHO/MDI images, 
Kuhn et al.l (120041 ) found no significant variations at any level above 7 mas, which is signifi- 
cantly lower than any variation reported from ground-based instruments (astrolabes) ranging 
from 50 to 2 00 ma s. Using the Solar Disk Sextant experiment aboard stratospheric balloons, 
Sofia et ajj (119941 ) indicated an antiphase variation between the solar radius and the 11-year 
cycle. 

Helioseismic methods can also be used to study the variations with the c ycle: in this 



case, a seismic rad ius is defined from oscillation frequencies of solar f-modes (ISchou et al. 



19971 ; lAntial Il998l ). The seismic radius is an interest ing probe of the subsu r face l ayers 
to a depth of about 15 Mm, as recently published by iDziembowski fc Goodei (120041 ) and 
Lefebvre fc Kosovichev (boosh (hereafter LK05). Using properties of /-mode frequencies, 
LK05 showed that the variations of the solar seismic radius were nonuniform with depth. 
More precisely, LK05 reported a non-monotonic change in the position of the subsurface 
layers: the radius changes between approximately 0.97 and 0.99 R & were in phase with the 
solar cycle while the radius of the shallower layers above 0.99 R & up to the surface changed 
in antiphase. In general, the results are consistent with previous conclusions t hat solar-cycle 
varia t ions in the solar radius are co nfined in the outermost layers of the Sun (lAntia fc Basu 
200i IDziembowski k Goodell2005h . 



Using an evolution code of the Sun, ISofia et al.l (120051 ) (hereafter SBDLT05) also found 
a non- homologous variation of the solar radius with depth and a radius change at the sur- 
face in a ntiphase with the so lar activity cycle. Using the Yale Rotating Evolution Code 
(YREC; (1 Winnick et al.ll2002l ) and including the effects of a variable magnetic field and field- 
modulated turbulence, SBDLT05 found a monotonic variation of the solar radius, sharply 
increasing in amplitude by a factor of approximately 1000 from a depth of 5 Mm to the solar 
surface between 1996 and 2000. 



1 the mas is the milliarcsecond and is equal to about 0.725 km on the Sun 
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Both, the helioseismology and model results show that the solar-cycle variations of the 
solar radius with depth in the subsurface layers are non-homologous, and that the radius 
change at the surface is in antiphase with the solar activity. 

However, there are also discrepancies between the model and helioseismic inferences 
concerning the amplitude of these variations and the behavior of the changes in the subsurface 
stratification in the first 15 Mm. LK05 found an oscillation and a change of phase in the first 
15 Mm with variations not exceeding 30 km, whereas SBDLT05 found a monotonic variation 
in the first 5 Mm with a factor approximately up to 1000 between the depth at 5 Mm and 
the surface. However, LK05 showed that the helioseismic inferences are not very sensitive 
to the surface changes. Therefore, to make a quantitative test of the model prediction in 
addition to the helioseismology data we use measurements of the solar surface radius from 
the ground and space. 



2. Computations and updated results concerning the inversion of /-mode 

frequencies 

As in LK05, data used in this study are frequencies of solar oscillation modes from 72 
day MDI observing rung§. Only /-modes for the period 1996-2005 have been selected. The 
oscillation modes, common for the whole period, are extracted to finally obtain 148 modes 
ranging from I = 125 to I = 285. We will here simply make a summary of the formalism 
used in LK05, and the reader is invited to see that paper for detailed explanations. 

A relation between the relative frequency variations 5v/ v for /-modes and the associated 
Lagrangian perturbation of t he radius 5r/r of subsurface layers has been established by 
Dziembowski fc Goode! J2004h : 



Su 



31 



2lu 2 I 



X 



where I is the degree of the /-modes, I is the moment of inertia, to is the eigenfrequency 
and g is the gravity acceleration. This equation is the starting point for our inversion that 
will permit us to obtain the variation of the position of t he subsurface layers determined by 
5r and Su/u. For these calculations, we used model S (jChristensen-Dalsgaard et al.lll996l ) 
calibrated to the seismic radius = 6.9599 x 10 5 km, and the standard Regularized Least- 



Square technique (jTikhonov Arseninl 119771 ) (as we have an ill-posed inverse problem). 



2 These files have been computed by J. Schou and are available on 
http : //quake . Stanford. edu/$\sim$schou/anavw72z/ 
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For the inversions, as in LK05, we select f- modes of angular de gree I below 250, for which 



near-surface magnetic and turbulence effects are not important (jLefebvre et al.ll2006l ) 



The inversion of Eq. [T] leads to the results presented in Figure [T] which is an update of 
Figure 3 of LK05. The figure presents the variation Sr for each year: every difference has 
been computed relative to the reference year of 1996. This updated graph presents the same 
features than the original one published in LK05: 1) No significant changes in the variation 
of the subsurface layers' depth below 0.97 R Q with a maximal amplitude of about 10 km at 
~ 0.985 R Q ; 2) Non-monotonic changes in the stratification, with the inner layer varying in 
phase with the solar activity cycle and the outer layer evolving almost in antiphase with the 
solar cycle with a maximal amplitude of about 26 ± 6 km at about 0.995 R Q . 

The new result is the radius changes are not precisely in antiphase with the 11-year 
cycle. The results for 2005 indicate there might be an additional phase lag between the 
radius change and activity. The 2005 data show the same behavior as 1997 data of the 
previous solar minimum but the new activity minimum is not reached yet. Significance 
of this result needs to be checked with future data; if confirmed it may have interesting 
implications. 



3. Comparison of model predictions with the data 

In this section, we compare our results with the model results obtained by SBDLT05. 
Figure 3 of SBDLT05 shows the radius change as a function of depth below the photosphere 
relative to the radius change at the depth of 5 Mm for their model 4, which includes both 
temporal variation of magnetic field and turbulence. However, the model does not provide 
an absolute value for the radius variations. Therefore, we have to calibrate the model to 
obtain the best match to the helioseismic f-mode measurements. 

The principle is to use Eq. [T]in a direct way by assuming a value of \AR 5 Mm\, obtaining 
— from model 4, and then calculating ^ from Eq. [TJ Then, we search for the | AR§Mm\ value, 
which provides the best match to the observe d changes of f-mode frequencies be tween 1996 



and 2000. For the reference we use model S (jChristensen-Dalsgaard et al.lll996l ) calibrated 



with the seismic radius (R@ = 6.9599 x 10 5 km), the same as in the YREC model. The 
first step was to digitize the curve at the top of Figure 3 of SBDLT05 and transform to the 
fractional radius. The panels of Figure [2] show: a) the digitized curve as given by SBDLT05; 
b) the AR computed with different \AR 5Mm \ ranging from 0.3 to 1 km; c) the associated 
AR/R; and finally, d) computed Avjv for the different AR/R (color curves) in comparison 
with the Avjv observed between 1996 and 2000 (points with error bars). 
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Only few sample curves are plotted in panel d) of Figure [21 and a study using the x 2 
parameter shows that the minimum of \ 2 is obtained for AR 5 Mm — 0.65 km, which gives 
a value for AR at the surface of approximately 600 km. However, such high variation of 
the solar radius has never been observed in direct solar limb observations. In particular, 



Kuhn et al.l (120041 ) put an upper limit of 7 mas, or just 5 km. This means that the theoret- 
ical model of the solar-cycle variations cannot satisfy both observational constraints, from 
helioseismology and from the limb measurements. To fit the f-mode helioseismology data it 
must have the surface ra dius variations by a factor 100 greater than the observational upper 
limit from SOHO/MDI (IKuhn et al.l 12004 ). In this case, it also exceeds the observational 
limit from ground-based measurements by a factor of 4-20. Vice versa, if the model radius 
is within the observational limit then the model cannot explain the f-mode data. 

Thus, we conclude that, according to the mathematical formalism previously described, 
the strong non-homologous variation of the subsurface solar radius proposed by SBDLT05 
is not compatible with the variations of the /-mode frequencies and the surface radius, 
observed between 1996 and 2000. These computations explain why our results published 
in LK05 and updated here are in contradiction in amplitude and behavior with the results 
of SBDLT05. Of course our computations are completely based on Eq. [TJ established by 
Dziembowski fc Goodd (120041 ) and on the hypothesis that the temporal variation of /-mode 
frequencies are mainly due to changes in the subsurface stratification during the solar activity 
cycle. On the other hand, the model proposed by SBDLT05 is based on the YREC code 
(I Winnick et al.l 120021 ) into which the effects of magnetic fields and turbulence have been 
included. But, as they said, the location, magnitude, and temporal behavior of the internal 
field are not known which imply to make assumptions on their treatment. We suggest that the 
difference obtained with our results could come from their treatment of the turbulence and 
magnetic fields in the subsurface layers. Overcoming this discrepancy between the theory 
and observations will help to understanding of the complicated magnetic and turbulence 
effects below the visible surface of the Sun. 



4. Discussion 



In this paper, we have shown that the solar-cycle variations of the solar radius that 
increase by a factor of approx imately 1000 from a depth of 5 Mm to the solar surface as 
published by ISofia et al.l (120051 ) (SBDLT05) are not consistent with the /-mode frequencies 
variations between 1996 and 2000 and the upper limit on the solar radius change from 
the simultaneous limb measurements. Using their model and the relationship between /- 
mode frequencies and solar radius variations, we obtained a variation at the solar surface 
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of about 600 km wh i ch is observed neither with ground-based nor space instruments. The 



model of ISofia et al.l (120051 ) includes the effects of a variable dynamo magnetic field and of 
a field- modulated turbulence with general assumptions on their magnitude and temporal 
behavior (based on an observed luminosity change of 0.1% and the shape of the activity 
cycle) because of the unknown location, magnitude and temporal behavior of these internal 
fields. Nevertheless, we are still cautious on our above results because (i) the near-surface 
effects of turbulence and magnetic fields are not treated in our approach and (ii) the lack of 
very-high /-mode degrees are necessary to understand the very outer layers (approximately 
above 3 Mm). So we do not exclude a bigger variation of the position of these layers than 
those plotted in Figure [1] and that the solar radius at the photosphere by limb observations 
(as measured by present instruments and future astrometric satellites) could be larger than 
the seismic radius issued from /-mode frequencies studies. Hence, we agree with SBDLT05 
and the necessity to measure precisely the variations of the solar radius from space. 
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x=r/Ro 



Fig. 1. — Radial variation (Sr) as a function of the fractional radius x = r/R®, obtained as 
a solution of the inversion of /-mode frequencies by a least-squares regularization technique. 
The reference year is 1996. The error bars are the standard deviation after average over a set 
of random noise added to the relative frequencies. The averaging kernels for this inversion 
are well localized between 0.985 and 0.996, with a typical half-width of about 0.003. 
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Fig. 2. — a): The ratio between radius change as a function of depth below the photosphere 
and radius ch ange at 5 Mm as a function of depth in Mm, i.e. same figure as in top panel 
of Figure 3 of Sofia et al. ( 2005I ) . b): Computed AR with different \AR 5Mm \ (see the legend 
expressed in km in panel c) as a function of the fractional radius x = R/R & . c): Computed 
AR/ R with different AR 5 Mm from legend (in km) as a function of the fractional radius x. d): 
Observed Avjv in black point with errorbars and computed integrated Avjv as a function 
of the degree I, using Eq. [Hand AR/R from left bottom panel. 



